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Abstract The MIP, a micro satellite probe onboard 
Chandrayaan -1 was designed to impact on moon surface at 
pre decided location for scientific data collection from close 
range which was released from Chandrayaan-1 mother ship 
on 14 th November 2008. This MIP impacted near the moon 
surface near south-pole after 24 minutes from its release and 
while descending MIP suffered a composite motions 
consisting of three fundamental motions viz. spinning, 
coning and forward translational (velocity). Such descending 
MIP showed a definite pattern of coverage on moon surface. 
Therefore in order to monitor and to select future landing site 
of impact on moon surface the associated component of 
motion viz.; spinning rate, coning rate, velocity of descend of 
MIP needs to be calculated and known. In view of this, an 
attempt has been made here to derive associated parameters 
of composite motions and attitude of MIP for predicting 
future impact site of descending probe on moon surface from 
the orbiting space craft using the time sequential images 
acquired by the descending MIP and the algorithm 
developed under MATLAB after Srivastava et. al. 2011. 
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I. Introduction 

Moon Impact Probe (MIP) , a micro satellite probe was 
the most important payload among the 11 payloads onboard 
Chandrayaan-1 designed by the Indian Space Research 
Organisation (ISRO) to impact the Shackleton Crator area 
near to south pole region of moon (Bussey et al 1999). The 
MIP took about 24 min to impact on the moon surface after 
it’s separation from the mother ship Chandrayaan-1 (Fig. 1) 
consisting of three other instruments viz.; Radar Altimeter, 
Mass Spectrometer and Video Imaging System ( Kumar 
2009,Goswami 2009) and while descending MIP has shown 


a definite pattern of coverage on moon surface. 

The mission sequence of the MIP is given as: 

T 0 -20 min MIP Power on 

T 0 -5 min MIP Enable Command 

T 0 MIP separation 

T 0 + 10 sec Spin up of MIP 

T 0 + 1489.02 sec MIP impact 

In present study the concept of fundamental motions of 
descending MIP has been applied for extracting individual 
motion parameters such as spinning rate , coning rate and 
velocity from the acquired sequential images of descending 
MIP over moon surface using the developed algorithm under 
MATLAB after Srivastava et al, 2011. 



Figure 1 . MIP Mission and path (Kumar ,2009) 


2. Data and Methodology Used 

In present study we have used near about 3000 - time 
sequential MIP photo images (720x576 pixels), acquired 
after the 25 minutes journey at 20:31 Hrs 1ST (Indian 
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Standard Time) covering the area around Shackleton Crater 
located in south polar region of moon and were received on 
14th of November 2008 by Chandrayaan-1. The time 
interval between two subsequent odd-even frames has been 
20 ms whereas it is 1.7 sec in between two subsequent 
even-odd frames ( Kumar 2009). These MIP images have 
shown a definite pattern of coverage on moon surface under 
the influence of composite motions of descending probe and 
from this pattern an attempt has been made to derive 
composite motions and the attitude parameters of MIP viz.; 
spinning rate, coning rate, velocity and inclination using 
algorithm developed after Srivastava et al., 2011. 
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Figure 2. Showing two subsequent odd and even frames. 


3. Analysis and Discussions 

A. Spinning Rate 

For calculating this rate two subsequent odd-even frames 
#VIS 2899 and #VIS 2900 ( Fig. 2) with time interval 20 ms 
have been taken where the point of rotation center has not 


changed and so by taking a common feature in both the 
frames one calculates the angle swept by the feature with 
respect to the center of rotation in a given interval of time. 
However the later frames of MIP (Fig. 3) have been taken 
when the view of footprints of MIP are almost vertical and so 
the distortion is minimum due to inclination. 

In above two frames the P and P’ features are same. In 
frame #VIS 2899, the pixel co-ordinate of P is (567,416) and 
the extended line from the rotation center cut the x -axis at 
point O (82,0). So the angle (O) subtended by point P with 
the x- axis is given by equation 1, 

® =tan _1 (416/567-82) =40.65° (1) 

Similarly ) in frame VIS # 2900 for point P’(532,453, the 
angle (O’) is given by the equation 2. O’ = tan" 1 

(453/532-147)= 49.65° (2) 

This means the features subtends an angle of 9° 
(=49.65°-40.65°) in 0.02 sec of time interval and in this 
fashion it completes one rotation in 0.8 sec ( = 0.02*360/9) 
which gives a spinning rate equivalent to 75 rotation per 
minute (RPM). 


B. Coning Rate 

The spinning motion of MIP causes in a definite pattern of 
coverage on moon surface where as the coning motion 
causes the features to reappear after some time in a set of 
images taken by the MIP. This could be seen clearly in initial 
frames where size of footprint is very large which reduces 
subsequently near the impact site ( Fig. 3) and thus 
possibility of revisiting of an object during one coning cycle 
is minimal. This is due to fact that the continued forward 
motion of camera displaces it away from earlier site before 
camera points back at the end of one coning cycle. 
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Figure 3. Variation in size of foot prints of MIP with inclination after 
separation. 

Reappearing objects have been observed and identified on 
a set of images and the time taken by that object to move 
from one extreme to the other and then return to its original 
position has been taken as one coning cycle period. Thus the 
frequency of coning cycle provides the coning rate of MIP 
motion. In the images shown below (Fig. 4) it is observed 
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that the features reappear / reorient themselves exactly after 
20 frames (ie 10 pairs of even and odd frames). 



Frame # VIS 1425 


Figure 4. Shows reappearance of objects exactly after 20 frames (#VIS 
1405 -#1425) 

In initial frame #1405, where the edge of landmass is 
almost horizontal which rotates itself by almost 180° after 10 
frames ie in frame # 1415. Finally the edge reorients as per 
the frame #1405 exactly after another 10 frames ie in frame # 
1425, however ,with some residual angel. In above three it 
also is observed that the ratio of deep space to land mass is 
changing from frame to frame and similar pattern is repeated 


after 20 frames. 

Now the time offset between 20 frames is 1.7*10 + 
0.02*10 = 17.2 sec which gives rotation speed equivalent to 
0.05814 rotation per second (RPS) and spinning rate 
equivalent to 3.485. This directly converts to a coning rate 
equal to 3.5 RPM (approximately) or its multiple which is 
the minimum possible value of coning rate. However it is not 
possible to measure accurately the rotation angle due to 
coning motion because it is possible to observe the effect of 
residual angle only at the end and hence it is not possible to 
know the number of rotation it has completed before it has 
reached to present position with residual angle and numbers 
of registered frames. Thus the minimum value of conning 
rate is 3.5 RPM approximately. 

Velocity-Forward Translational Motion 

For this we have taken two frames (# VIS 2871 and #VIS 
2874) where the object (P) has appeared and the view of the 
foot prints of MIP are nearly vertical with angle is less then 
4° as shown in Fig.5. 


Table 1. Details of Frames 


Frame # 

View 

angle 

(degree) 

Altitude 

(km) 

Foot print 
(Km 2) 

Pixel no. 

2871 

2.3297 

21.18 

10.93x8.44 

720x576 

2874 

2.2776 

20.97 

10.82x8.36 

720x576 


Here the total time difference (t) between these two frames is 
1.74 secs as explained below. 

#2871— 20msec — # 2872— 1 ' 7sec — # 2873~ 20msec — #2874 
So the elapsed time in between the above two frames is 
1.74 sec and by taking the calculated coning rate as 3.5 RPM 
(=0.05833 rotation per second ) and the spinning rate as 72 
RPM (= ie 1.2 rotation per second ) the total rotation comes 
out to be 36.42° and the total rotation angle due to spin equals 
to 751.68°(=1.2* 1.74). Thus the total rotation due to coning 
and spinning comes out to be 88.1°(=36.42° + 751.68° ) 
which could be attributed to two full rotation (=720°) plus a 
residual angle of 68.1°. 



Frame # VIS 2871 
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Frame # VIS 2874 


Figure 5. P 0 and P f shows the position of objects in frame #VIS 2871 and # 
VIS 2874 respectively with almost vertical view due to coning and spinning 
motions of MIP. 

In frame # VIS 2871, the co-ordinate of the feature is 
P 0 (x 0 ,y 0 ) =("208,53) and in the frame # VIS 2874 the 
co-ordinate of final position is P f (x f ,y f ) = (66,-153) and such 
displacement in position of features has happened due to 
composite motions of MIP viz. the spinning, coning, and 
translational. 

Therefore by applying the rotational matrix R with 
residual (0) = 68.1° on the initial co-ordinate P 0 , one get the 
amount of displacement due to the spinning and coning only. 
The difference between the final position, P f and the 
position after rotation shall give the effect of the velocity 
component only for which the pixel coordinate has comes 
out to be ,P S (-127.05,-172.95). 

This point P s represents the position of initial point, P 0 
(x 0 ,yo) under the influence of two motions of MIP viz.; 
spinning and coning without the effect of velocity. Therefore 
the difference in position in between the point P s and the 
final point P f shall give translational component of initial 
point features due to the velocity. This distance (d) between 
these two points P s and P f is approximated to be around 194 
pixel as derived as shown below using equation (3) . 

d = V((x f .x) 2 +(y f _yf) -(3) 

= 194.1 pixels 

and the rotational angle between these two points is given by 
equation (3). 

a = tan" 1 ( y f _y )/(x f . x’)= 5.92 0 = 6° approx (4) 

The distance, d, is the distance between the final point P f 
and the point P s after rotation due to the spinning and coning. 
However we know that the MIP has travelled a distance in a 
time interval, a translational (forward) motion giving the 
velocity of motion. Hence the distance travelled by the object 
in the frames is only due to velocity and the calculation of 
which is described below: 

Taking the reported velocity value of MIP as 1.6363 


km/sec (Kumar 2009), the distance travelled by the MIP 
comes out to be 2.8471 km. in time interval of 1.74 sec and 
then considering the average foot print size (720x576) pixels 
the corresponding size of one pixel equals to 0.015 km (= 
10.85 km/720) on moon surface. Thus in fact the travelled 
distance (d) in terms of pixel dimension (p) is 189.8 pixels. 

Thus from above calculation we observed that the camera 
in a given time has moved from point P st 0 P’ and is 189.8 
pixel away. However the calculated distance in terms of 
pixel between the shifted position, P s and the final position , 
P f has come out to be 194.1 pixels with orientation of motion 
vector as 6°. But after de rotating the images (Fig. 6) the 
angel of orientation of motion vector comes out be 57°. 
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Figure 6. The space sequential images and direction of motion vector after 
de-rotation. 

Now we are observing that the direction of velocity vector 
is not matching with the value derived from present study. 
Probably we have considered small value of coning rate 
which may be multiple of cycle and so the sampling rate is 
very coarse which is not adequate to resolve the accurate 
number of rotation (n). 

4. Conclusions 

Using the time sequential images acquired by the 
descending probe and the resulting pattern of surface cover 
on moon it has been possible to estimate parameters of 
motion components associated with descending MIP viz.; 
the spinning motion, coning motion and velocity (forward 
translation motion) which are described below: 

1. From the sequential images, the spinning rate of the 
MIP has been obtained as 75 RPM where as the system 
Accelerator has provided a value of 72 RPM . 

2. The derived value of conning rate of MIP is 3.5 and its 
multiples. It has become difficult to determine distinct 
number of multiples as because the numbers of registered 
frames are not adequate to calculate this. 

3. The distance travelled by the MIP is 2.9115 km 
(194.1*0.015) which gives a velocity equivalent to 1.6732 
km/sec as against the reported 1.6363 km/sec. velocity due to 
reason as explained above.. 
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Such discrepancies in calculated motion parameters are 
probably occurring due to association of an external force 
which modifies the motion of MIP and nature of which is to 
be ascertained. 


Acknowledgments 

Authors express their sincere thanks to the Director, Space 
Application Centre (ISRO) Ahmedabad, India and also to 
Director, Indian School of Mines, Dhanbad, India for 
providing required facilities to carry out above work. Junior 
author (TKN) in particular is thankful to Scientists of Space 
Application Centre, Ahmedabad for permitting him to work 
as Intern on this problem. 


REFERENCES 

[1] Bussey DBJ, BD Spudis and MS Robinson (1999), 
“Illumination condition at the Lunar South pole”, Geophysics 
Research Letter ,26, pp 1187-1190. 

[2] Goswami J.N., M. Annadurai, (2009), “Chandrayaan-1: 
India’s first planetary science mission to the moon,” Current 
Science, vol. 96, no.4. pp.486-491 

[3] Kumar Ashok and MIP Project Team, ( 2009), “The Moon 
Impact Probe on Chandrayan -1”, Current Science, Vol 96,No 
4 Feb 2009 pp 540-543. 

[4] Srivastava V K, Nandy T K and Verma A, (2011), “Study of 
motion components associated with descending Moon 
Impact Probe (MIP) to moon surface from the chandrayaan-1 
space craft from simulated surface motion traces using the 
developed algorithms under MATLAB and from sequential 
photo images acquired during the descend”, Int. Journal of 
Experimental Algorithm ,Vol 2, Issue 3,pp 74-81. 



